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With an eye toward LHC processes in which theoretical precisions of 1% are desired, 
we introduce the theory of the simultaneous YFS resummation of QED and QCD to 
compute the size of the expected resummed soft radiative threshold effects in precision 
studies of heavy particle production at the LHC. Our results show that both QED and 
QCD soft threshold effects must be controlled to be on the conservative side to achieve 
such precision goals. 
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1. Introduction 

At the LHC/ILC, the precision requirements for soft multiple gluon (n(g)) effects 
will be even more demanding than at FNAL, where the uncertainty on mt [1], 
5m t = 4.3 GcV, receives a soft n(g) uncertainty ~ 2-3 GeV, and soft n(g) MC 
exponentiation results will be an important part of the necessary theory - YFS 
exponentiated 0(a 2 s )L calculations, in the presence of partem showers, on an event- 
by-event basis. 

As many authors [2] prepare the necessary results that would lead to such a 
precision on QCD for LHC processes, the QED corrections need to be addressed as 
well. Estimates by Refs. [3-7] show that one gets few per millc effects from QED 
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corrections to structure function evolution. In this paper, estimate the size of QED 
corrections at threshold at the LHC. 

Treating simultaneously QED and QCD in the respective YFS [8,9] exponenti- 
ation we discuss threshold effects at the LHC in the candidate luminometry [10,11] 
processes pp — > V+n(j)+ni(g)+X — > it' +n'(j)+m(g)+X , where V = W , Z,a.nd 
I = e,fi, £' = u e , Vn(e, fx) respectively for V = W + (Z), and I — v^v^, t' = e, jx 
respectively for V = W~ . 



2. YFS Theory and its Extension to QCD 

In Refs. [9] the renormalization group improved YFS theory [12] for e + (pi)e~(qi) — > 
f(P2)f(q2)+n(j)(k 1 , •, fc„) is realized by Monte Carlo methods, where the respective 
cross section da exp and all of the attendant IR functions, {B, B, D, S}, and hard 
photon residuals ,{/3„(fci, . . . , fc„)}, are specified in Refs. [9]. In Refs. [13,14] we 
have extended the YFS theory to QCD: the net result is that in the analogous YFS 
theory we have the replacements 2aReB + 2a B — > SUMir(QCD), D — > Dqcd, 
and f3 n (ki, . . . , k n ) — > j3 n (ki, . . . , k n ), where the QCD YFS functions are defined in 
Rcf. [13] and the gluon residuals [13] /3„(fci, . . . ,k n ) are free of all infrared diver- 
gences to all orders in a s (Q). The genuine non-Abclian IR physics is encoded [13] 
here in the j3j . The YFS resummation which we discuss here is fully consistent with 
that of Refs. [15,16]. See Ref. [17] for more discussion of this point. 



3. Extension to QED&QCD and QCED 

Simultaneous exponentiation of QED and QCD higher order effects gives [17] 



da cxp = gSUMirt(QCED) £ / TJ ?_5Z! ±±21 
m,n— j± — 1 J2=l 



d 3 k 3l fr fk^ r _d*y_ 
*Ji M k' h J (2^)4 



p iy(Pi+9i-P2-«2-E feji-E fe' J2 )+- D QCED / 5 /r. . ./ / . d 3 p 2 d 3 q 2 

P 2 32 

(1) 

where the new YFS functions, 

SUMui(QCED), -Dqced and j3 m:n (ki, . . . , k m ; k[, . . . , k' n ), where the latter has 
m hard gluons and n hard photons, are defined in Ref. [17]. The infrared algebra 
QCED [17] obtains: the average Bjorken x values for the QED and QCD emissions 
imply [17] that QCD dominant corrections happen an order of magnitude earlier 
than those for QED so that the leading -level gives a good estimate of the size 
of the effects we study. 

4. QED®QCD Threshold Corrections at the LHC 

For the basic formula (we use the standard notation here [17]) 

da exp (pp -> V + X — ► U! + X') = ^ J dx l dx J F l (x l )Fj(x : j)daexp(x l x : js), (2) 

i,3 
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we use the result in Q for V — Z here with semi-analytical methods and structure 
functions from Ref. [24]. See also the work of Refs. [18-23]. A Monte Carlo real- 
ization will appear elsewhere [25], wherein we will ultimately use HERWIG [26], 
PYTHIA [27] and /or the new shower algorithm in Ref. [29] in lieu of the {F,} 
and thereby, in principle, improve on the shower/exact result combination in 
Ref. [28] . Due to its lack of the appropriate color coherence [30] , we do not consider 
ISAJET [31] here. 

We compute , with and without QED, the ratio r exp = o~ exp / 1 <JBom to get the 
results (We stress that we do not use the narrow resonance approximation here.) 



QED is at the level of .3% at both LHC and FNAL a . We agree with the results 
in Refs. [18-22]. The QED effect is similar in size to structure function results in 
Refs. [3-7]. 
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